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A sample of over 1000 objects selected from a 1.4 GHz survey made by 
the Australia Telescope Compact Array (ATCA) is used to study the proper- 
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""^ 1 deg 2 , is 50% complete to 0.2 mJy. Over 50% of the radio sources are found to 

have optical counterparts brighter then R rj 21.5. Spectroscopic observations 
2 ! °f 249 optically identified radio sources have been made, using the 2-degree 

Field (2dF) facility at the Anglo- Australian Telescope (AAT). Redshifts and 
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equivalent widths of several spectral features (e.g., Ha and [OIIJA3727) sen- 
sitive to star formation have been measured. On the basis of the photometric 
and spectroscopic data, the optically identified radio sources are classified 
as (i) absorption-line galaxies, (ii) star-forming galaxies and (iii) Seyfert-like 
galaxies. 

The spectroscopic sample is corrected for incompleteness and used to esti- 
mate the 1.4 GHz and Ha luminosity functions (LFs) and luminosity density 
distributions. The 1.4 GHz LF of the star-forming population has a much 
steeper faint-end slope (1.85) than the ellipticals (1.35). This implies an in- 
creasing preponderance of star-forming galaxies among the optically identified 
(i.e., z <1) radio sources at fainter flux densities. The Ha LF of the faint radio 
population agrees with published Ha LFs derived from local samples selected 
by Ha emission. This suggests that the star-forming faint radio population is 
coincident with the Ha selected population. 

The 1.4 GHz and Ha luminosity densities have been used to estimate the 
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star formation rates (SFRs). The two SFRs agree, both giving a SFR density 
of 0.032 M yr _1 Mpc~ 3 in the range z <1. Radio selection appears to be as 
effective as Ha selection in finding the galaxies which dominate star formation 
at a given epoch. Although the sample contains many galaxies lying beyond 
z w 0.3, it does not reveal a significant rise in the global star formation rate 
with increasing redshift. This result suggests that the optical counterparts of 
galaxies undergoing vigorous star formation at redshifts beyond z 0.3 are 
generally fainter than R w 21. 

Key words: cosmology: observations - galaxies: starburst - galaxies: star- 
burst - radio continuum: galaxies 



1 INTRODUCTION 

The population of apparently bright radio sources (Si. 4 > 100 mJy) is dominated by classi- 
cal radio galaxies, whose optical hosts, when detectable, are usually identified as luminous 
elliptical galaxies with red colours and absorption-line optical spectra, although some are 
hosted by QSOs. The elliptical hosts sometimes emit nebular emission lines, which can be 
quite strong for the most luminous radio sources ( [Hine fc Longair 1979| ). The lines appear 
to be excited by an active galactic nucleus (AGN) radiating a non-stellar continuum ( [Baum 
fc Hcckman 19891) . The normalised source count distribution of the classical radio galaxies 
falls sharply with decreasing flux density, as a result of cosmological dimming and strong 
source evolution for z < 2 ( |Wall fc Jackson 19971 ). 

The apparently bright radio source population contains, in addition to these classical 
radio galaxies, a small proportion of spiral/irregular galaxies. These are markedly different 
from the classical radio galaxies. For example, even samples selected to have large radio-to- 
optical luminosity ratios have a radio luminosity smaller than that of most classical radio 
galaxies QCondon et al. 1982|) . The core/jet/lobe radio morphology of classical radio galaxies 
is absent, and replaced by centrally concentrated but nevertheless diffuse emission. More- 
over, a high proportion of radio-bright spirals have peculiar optical morphologies and strong 
emission lines excited by stellar photospheric radiation, apparently reflecting an enhanced 
rate of star formation flCondon et al. 1982| ). 

Spiral/irregular galaxies do not contribute in significant numbers to radio source counts 
for Sxa ^ 10 mJy. However, radio surveys reaching fainter flux density limits, S1.4 < 1 mJy, 
reveal a significant change in the slope of the normalised source count distribution at 
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S 1A ps 5 mJy [Windhort, van Heerde & Katgert ( |1984| ); Condon ( |1984| ); Fomalont et al. 
(|98D; Mitchell & Condon (p85|) ; Windhorst et al. fljgSg )]. The change suggests the possi- 



ble existence of a numerous population of apparently faint radio sources unlike the classical 
radio galaxies. Optical photometry and spectroscopy of the population reveals that many 
of the faint radio sources are associated with blue galaxies often exhibiting peculiar (com- 
pact, interacting and merging) morphologies. Kron, Koo & Windhorst (|1985|) and Thuan & 
Condon ( |1987|) have shown that these galaxies are quite unlike the hosts of classical radio 



galaxies, but are reminiscent of radio-bright local spiral galaxies. 

The central astrophysical importance of the faint radio population lies in the fact that the 
unevolving, local luminosity function of the non-AGN radio source population can explain 
neither the large number counts nor the source-count distribution at sub-mJy flux densities 
( Condon 1984|) . Attempts to resolve this problem show that significant evolution may have 



occurred over the redshift range spanned by the observed population. For example, Danese et 
al. ( P-987| ) established consistency by adopting a model which includes luminosity evolution 
of a starburst /interacting population over an e-folding time of about 25% of the Hubble time. 
There are hints that the IRAS galaxy population undergoes similar evolution, although the 
IRAS surveys did not extend to the faint levels corresponding to deep radio surveys ( |Lonsdali 



et al. 19871) . Population evolution of FIR-selected galaxies has been confirmed by deeper FIR 



observations made by the Infrared Space Observatory (ISO) flRowan-Robinson et al. 1997 ). 



Additional evidence, consistent with the hypothesis of recent (z < 1) rapid evolution of a 
large fraction of the faint radio population, has come from the improved redshift distribution 



statistics for the optical counterparts of faint radio sources reported by Benn et al. (|1993| ) 
and Hopkins et al. ( |1999j ). 

The cosmological significance of this evolving population has been clarified by opti- 
cal/UV studies using deep ground-based observations and data from the Hubble Space Tele- 
scope (HST). Studies by Madau et al. (|1996|) imply not only that the space density of 
star- forming, morphologically disturbed galaxies increases with redshift to z > 1, but also 
that this increase is associated with a corresponding increase in the global star formation 
rate in the Universe. A similar picture emerges from observations of faint sub-mm sources 
( |Blain et al. 1998| ). 

The 1.4 GHz luminosity of a spiral/irregular galaxy can be used to estimate its star 
formation rate (SFR), assuming the appropriate calibration (Condon & Yin 1990; Condon 
1992). Using this result and published 1.4 GHz luminosity functions, Cram (|1998|) has argued 
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that the global rate of star formation appears to be dominated by galaxies with SFR > 15 
M yr" 1 . This implies that galaxies with Pi. 4 > 5 x 10 22 W Hz" 1 are likely to dominate 
the global SFR. Such a galaxy would be brighter than 10 /zJy at z — 1, implying that deep 
radio surveys probe the population which dominates global star formation to cosmologically 
significant redshifts. 

The discoveries summarised above suggest that the faint radio source population may 
well contain all of the galaxies which dominate the global star-formation processes revealed 
by optical/UV diagnostics, at least for redshifts out to z ~ 1 at the current limits of radio 
surveys. However, the implications of this have been explored directly for only a few dozen 
objects [Benn et al. (|1993|) ; Hammer et al. (|1995|) ]. In the present study, we use a sample of 
over 420 faint, radio-selected galaxies with optical data to define a well characterised sub- 
sample of star-forming objects, many (over 50%) with measured redshifts. From this sub- 
sample, we use the radio and Ha luminosities to estimate luminosity functions, luminosity 
densities and star- formation rates in the era z < 1. We show that the radio and optical 
properties of this sample are consistent with one another, and with the properties of samples 
from the local Universe. We use the results to provide an unbiased estimate of the star- 
formation history (free from dust effects) in the surveyed volume. 

In the following section we summarise the observations. Section 3 presents the luminosity 
functions and star formation rates, and Section 4 discusses the results from this investigation. 

2 OBSERVATIONS 

The observations have been accumulated during several campaigns connected with the 
Phoenix Deep Survey. The criteria used to select the field, as well as the radio observations 
and data reduction methods, are discussed by Hopkins et al. (|1998|) . The radio observations 
were made with the Australia Telescope Compact Array (ATCA) at 1.4 GHz, covering an 
area of diameter 2°, centered at RA(2000) = 01 h 14 m 1216 ; Dec. (2000) = -45° 44' 8"0. A 
total of 1079 radio sources are detected to S1.4 ~ 0.1 mJy The radio flux distribution for the 
entire survey is presented in Fig. 1 (solid line). The survey is 60% complete to Si. 4 ~ 0.3 mJy, 
with a total of 804 sources detected to this flux limit. The completeness drops to 50% at 
S1.4 ~ 0.2 mJy, with a total of 964 galaxies to this limit. The 80% completeness limit is 
achieved at about S1.4 = 0.4 mJy with 656 galaxies detected. In the present study, we con- 
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Figure 1. 1.4 GHz radio flux distribution of detected radio sources (solid line), radio sources with optical identifications (dotted 
line) and radio sources with spectroscopic measurements (dashed line). 

sider galaxies with Si. 4 > 0.2 mJy and refer to this in the following sections as the 'complete' 
radio sample. 

CCD observations of the entire radio survey area were attempted in the V and it- 
bands, using the Anglo- Australian Telescope (AAT). The observations were performed to a 
magnitude limit of R = 22.5, corresponding to the limit expected for successful spectroscopy 
with the two-degree field facility (2dF) of the AAT. Over 50% of the radio sources have 
optical identifications, corresponding to a total of 420 galaxies to S1.4 ~ 0.3 mJy and 490 
galaxies to 0.2 mJy. The radio flux distribution for galaxies with optical IDs is shown in Fig. 
1 (dotted line). The optical photometry and radio-optical source identifications are described 
in Georgakakis et al. (|1998|) . 

Spectroscopic observations were carried out for a random sample of the optically iden- 
tified radio galaxies to a limiting magnitude of R = 21.5, using the 2dF facility. Accurate 
redshifts were successfully measured for a total of 228 radio sources, using several lines, 
including H/3, Ha, [OII]A3727 and [OIH]A5003 flGeogakakis et al. 1998Q . Of this number, 
210 galaxies are in the 'complete' radio survey (i.e., have Si. 4 > 0.2 mJy), implying « 50% 
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completeness for the spectroscopic sample used here. The fraction of galaxies with spectro- 
scopic data at any given radio flux is presented in Fig. 1 (dashed line). Diagnostic spectral 
line ratios were estimated from these spectra and used to classify the galaxies in the spectro- 
scopic sample. This is explained in detail in Georgakakis et al. (1998) where the galaxies are 
classified as (i) elliptical, (ii) starforming, (iii) Seyfert and (iv) a small number of objects for 
which the spectral quality was not high enough to allow an accurate classification. Examples 
of the 2dF spectra of the radio sources classified to the above four sub-types are shown in 
Fig. 2. 

Spectroscopic signatures of probable star formation (mainly [OII]A3727 and Ha) are 
detected in the spectra of 63% of the faint radio sources. The redshift distributions for the 
complete radio sample (Si. 4 > 0.2 mJy) is compared in Fig. 3 with that of the sub-sample 
showing evidence of star formation activities as revealed from their spectra. The majority 
of star-forming radio sources brighter than the optical magnitude limit of this survey have 
^<0.3. 

One aim of this study is to derive quantitative optical measures of the star formation 
rate, using Ha wherever possible, to compare with the results from the radio data. However, 
an absolute flux calibration is not available for our 2dF spectra. Thus, the Ha luminosity 
(Pho) °f eac h galaxy has been estimated from its measured Ha equivalent width, EW a , and 
the i?-band magnitude. The spectral synthesis models of Bruzual and Chariot (|1993|) were 
then used to estimate the .R-band i^-corrections and the rest-frame continuum luminosity 
at Ha, C a ( Gcogakakis et al. 1998 ). The Ha luminosity is then estimated using 

P Ha = C a x EW Ha . 

For the few galaxies with z > 0.3, Ha shifts out of the observed spectral range. For these, 
we used [OII]A3727 as a proxy for Ha, converting its equivalent width using the relation 
( IKennicutt 1991 ) 

EWon = 0.4 x EW Ha . 

As noted by Kennicutt (|1992| ), the scatter about this correlation has an RMS dispersion of 
about 50%, due in part to galaxy-to-galaxy differences in extinction. 
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Figure 2. 2dF spectra for the 4 types of galaxies classified here. The identified spectral lines arc marked with the estimated 
rcdshift of the objects specified. The spectra correspond to; Elliptical (bottom left), Star-forming (top left), Seyfert 2 (bottom 
right) and 'unclassified' (top right) 

3 LUMINOSITY FUNCTIONS AND STAR FORMATION RATES 

In this section we estimate the 1.4 GHz and Ha luminosity functions for our sample. These 
will then be used to determine the luminosity density distribution and, for the star-forming 
galaxies, the star formation rate. To ensure an adequate sample size, we first assume that the 
sample statistics are independent of redshift (i.e. that there is no dependence of the selection 
functions on redshift). We consider two radio samples, consisting of the elliptical and the 
star-forming populations, as classified from their optical spectra (the small proportion of 
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Figure 3. The redshift distribution (solid line) of the Phoenix sample compared with the distribution of galaxies with Ha (or 
[OII]A3727) detections (dotted line). 

Seyfert or unclassified galaxies is excluded). Our Ha sample combines all galaxies in the 
radio-selected sample which exhibit Ha emission, since this would appear to provide a sample 
selected in much the same way as the objective prism sample of Gallego et al. ( |1995| ). 

3.1 Sample completeness and selection effects 

There are two sources of incompleteness in the present spectroscopic survey. Firstly, only 
radio sources with optical counterparts have been observed spectroscopically. This biases the 
spectroscopic sample against optically faint galaxies (e.g., distant, dusty and/or low surface 
brightness objects). Secondly, only a randomly selected sub-sample of the optically identified 
galaxies which satisfy the optical magnitude range accessible by the 2dF have been observed 
spectroscopically. 

To correct for the first bias we define a completeness function at any given radio flux 
density, T]i (S1.4), as the ratio of the number of radio sources with optical identifications 
to the total number of sources in the complete radio survey. This completeness function 
is displayed in Fig. 4a. The second bias can be corrected using a completeness function at 
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Figure 4. (a) The selection function r}\ (Si. 4) defined as the ratio of the number of faint radio sources with optical identifications 
to the total number of radio sources in the complete radio survey, (b) The selection function r)2(K) defined as the ratio of the 
number of galaxies with measured spectra to the total number of sources with optical counterparts in the complete radio survey. 

any given optical magnitude, r) 2 (R), defined as the ratio of the number of galaxies with 
measured spectra, to the total number of sources with optical counterparts in the 'complete' 
radio survey (S1.4 > 0.2 mJy). This completeness function is shown in Fig. 4b. The final 
completeness function is then pi(Si A ,R) = 771(61.4) ^(R). 

The Ha sample systematically excludes optically identified but faint galaxies which yield 
spectra with low S/N ratios and hence no reliable detection of Ha (or [OII]A3727). A com- 
pleteness function to correct for this selection effect, p 2 (R), will depend on the optical mag- 
nitude of an identification, and is defined as the ratio of the number of objects with Ha 
detection to the total number of sources with optical identifications in the complete radio 
survey. The function p 2 (R) is shown in Fig. 4c. 
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Figure 4. (c) The selection function of the Ha survey, p2(R), defined as the ratio of the number of galaxies with Ha detected 
to the total number of galaxies with optical identification. 



3.2 The 1.4 GHz and Ha luminosity functions 

Using the completeness functions derived above, we correct the sample for incompleteness 
and selection effects and estimate its radio [$1.4(^1.4)] and Ha [$ a (Pffa)] luminosity functions 

as 



where the summation is over the number of galaxies in radio or Ha flux bins in the complete 
sample, with V max given by 



Here, is the completeness function corresponding to the radio (k = 1) and H Q (k = 2) 
surveys, is the luminosity distance, and the integral limits are z 1 = max(z ra d, z opt , z min ) 
and Z2 = min(z' rad , z' opt , z max ) ■ The quantities z raa > and z opt are, respectively, the redshifts 
corresponding to the bright radio flux and optical i?-band magnitude limits of the survey 
while z' rad and z' opt are the maximum redshifts at which a source could be included (at the 
flux limit of the survey) in the radio sample and be detected at the optical wavelength. The 
quantities z m i n and z max are respectively the minimum and maximum redshifts covered by 
galaxies in the radio survey, i^-corrections are estimated in the i?-band using model spectral 
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energy distributions for intermediate-type spirals (section 2), and at 1.4 GHz assuming an 
spectral energy distribution proportional to z/ -0 7 . 

A parametric form which successfully models the luminosity functions of the IRAS ( |Saun- 
ders et al. 1990| ) and faint radio ( [Rowan- Robinson et al. 19"93| ) sources has been adopted to 
fit the luminosity functions in this study. It is the function 

$(P) dilog P) = C*(^y- a exp[- JL/o^ o( i + 21)] d{ i og P)j 

where C*, P*, a and a are the fitted parameters. Fig. 5 displays the radio (1.4 GHz) LF for 
all the galaxies in the 'complete' radio survey (Fig. 5a), the star-forming population (Fig. 
5b) and the ellipticals (Fig. 5c). The error bars represent uncertainties assuming Poisson 
statistics. The < V/V max > estimates are given in Table 1. The parametric fits to the radio 
LFs, assuming the above form, are carried out for cases with no evolution and assuming 
luminosity evolution in the form Pi A (z) = Pi.4(0) (1 + z) 3 . This form of luminosity evolution 
is adopted to be consistent with recent optical (Lilly et al 1996), far- infrared (Kawara et al 
1998) and radio (Rowan-Robinson et al 1993) observations. The fits were performed to the 
LFs corresponding to all the galaxies in the 'complete' radio survey, the starburst population 
and elliptical galaxies. Figure 5 compares the results for the two cases, assuming no-evolution 
(solid line) and luminosity evolution (dotted line). The estimated 1.4 GHz LF parameters 
are listed in Table 1. For a given population of galaxies in Figures 5a, 5b and 5c, we find 
no significant difference in the radio LFs between the two cases (assuming no evolution 
or luminosity evolution in the form (1 + z) 3 ). This result argues against the presence of 
significant luminosity evolution in the radio LFs in the redshift range, < z < 1, covered 
by the present survey. 



The LF for the star-forming radio sources in Fig. 5b is compared with Condon's ( 1989 ) 
estimates from the bright spirals (open circles), starbursts (filled triangles) and a sample 
of UGC galaxies with diameter > 1' (asterisks). The present sample has a mean redshift 
of < z >= 0.3, compared to < z >= 0.05 for galaxies from Condon (1989) and Sadler et al 
(1989). Therefore, the comparison between these samples in Figures 5b and 5c assumes no 
evolution for the radio sources in the range < z < 0.3, as justified from the above results. 
The independently derived LFs show excellent agreement in both shape and normalisation. 
Moreover, extrapolation of our radio LF in Fig. 5b to faint fluxes (Pi. 4 < 10 21 W/Hz) agrees 
well with the faint-end of the radio (1.4 GHz) LF for the spiral and irregular galaxies in 
Condon ( |1989| ) . Fig. 5c compares the radio LF derived here with that of the elliptical galaxy 
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Figure 5. (a) The 1.4 GHz luminosity functions for all galaxies in the 'complete' radio survey (filled circles). The sample is 
corrected for incompleteness. The curve is a parametric fit to the LF, using the function described in the text and assuming 
no evolution (solid line) and luminosity evolution in the form (1 + z) 3 (dotted line). . The LF is estimated in the range z < 1. 



radio luminosity function from Sadler, Jenkins & Kotanyi ( 1989 ) converted to 1.4 GHz 
assuming S u oc v . Considering that these samples have different selection criteria, depths 
and completeness limits, the agreement between the two LFs in Fig. 5c is remarkably good. 

The Ha LF of the sub-mJy radio sources with Ha detection is shown in Figure 6. The 
error bars are estimated assuming Poisson statistics. Fits to the analytical form discussed 
above are also shown, and the resulting parameters are listed in Table 2. The Ha LF agrees 
with that determined from an Ha survey of field galaxies ( Gallego et al. 199"5| ). (The faintest 
point, at log(Pff a ) = 33.5 W, significantly deviates from the fitted LF in both studies and 
is therefore excluded). Applying a Schechter LF fit to the present sample gives essentially 



the same parameters as those derived in (|Gallego et al. 1995|) , implying that the result 
here is independent from the parametric form of the LF. Considering the selection effects 
in the present survey (Section 2) the agreement is reassuring, and implies that Ha-emitting 
objects in the faint radio population are similar to the normal Ha-emitting galaxies in the 
local Universe. 



© 0000 RAS, MNRAS 000, 000-000 



Star formation in faint radio sources 13 




22 23 24 

log (P 14 W Hz- 1 ) 



26 





-2 




-3 


X 




CD 




-d 






-4 


u 




ft 




id 






-5 






o 





_ Elliptical Galaxies 


(c) - 


■ ' ' f _ f 


\. 4 - 



20 21 22 23 24 

log (P I4 W Hz-' ) 



25 



Figure 5. The 1.4 GHz luminosity functions for the faint radio population, corrected for incompleteness. Error bars are 
estimated assuming Poisson distribution. The curves present parametric fits to the LFs, using the function described in the 
text and assuming no evolution (solid lines) and luminosity evolution in the form (1 + z) 3 (dotted lines). The LFs are estimated 
in the range z < 1 and correspond to (b) star-forming galaxies here (filled circles) compared to the spiral/star-forming (open 
circles) and starburst (filled triangles) galaxies from Condon (1984) and (c) elliptical galaxies from this study (solid circles) 
compared with ellipticals from Sadler et al. (1990)-(open circles). 

The slope of the faint-end of the radio LF increases from a = 1.35 for ellipticals (Fig. 5c) 
to a = 1.85 for the star-forming population (Fig. 5b). A similar result was also reported by 
Condon ( 1989|) . It implies that in the redshift range < z < 1 and at modest radio powers 
(Pi.4 < 10 23 W/Hz), there is an abundance of star-forming galaxies compared to ellipticals. 
At the median redshift of our survey, z « 0.3, galaxies with P L 4 w 10 23 W Hz -1 have 
S1.4 ~ 0.2 mJy, and most of them have optical counterparts which will have been detected 
optically at a limit of R ~ 22. Thus, we may conclude that for the optically identified part 
of the sub-mJy population, the number of star-forming "IRAS-type" galaxies exceeds the 
number of ellipticals. Further observations of samples which are complete to known radio 
and optical limits would be of great value in elucidating the connection between the shape of 
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Table 1. Parameters of the analytical fits to the 1.4 GHz luminosity functions and the < V/Vmax > estimates. The results arc 
presented for all the faint radio sources in the present survey and for the starforming and elliptical sub-classes. Two different 
cases, assuming no evolution and luminosity evolution in the form Pi. 4(2) = Pi.4(0)(l + z) 3 arc considered. 
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Table 2. Parameters of the analytical fits to the Ha luminosity functions and the < V/Vmax > estimates. The results are 
presented for all the faint radio sources in the present survey and for the starforming population 
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the LF for z < 1, possible evolutionary effects, and the observed source count distributions 
including objects which have not yet been optically identified. 

There are seven detections of elliptical galaxies with Ha in the present survey. It is 
important to note that sub-mJy radio sources classified as elliptical galaxies have rest frame 
V — R colours close to those of 'normal' ellipticals (V — R w 0.6) and follow their V — R 
colour-magnitude relation (Georgakakis et al. 1998). Therefore, although the radio emission 
may be from an AGN, there is little evidence of AGN emission in the optical continuum. 
It is not excluded that both the radio and the Ha emission from these objects is related at 
least in part to star-forming activity at a level too weak to dominate the optical continuum 
colours. 

3.3 Luminosity density distributions and star formation rates 

The main source of uncertainty in converting the observed luminosity to star formation 
rate (SFR) is the dependence of the calibration on a number of poorly known parameters, 
including the initial mass function and the history of star formation prior to the epoch 
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Figure 6. The Ha luminosity function with a parametric fit (solid lines) for all galaxies with Ha detectio n in the radio surve y. 
This is compared with the Ha luminosity function from an H a-selected sample of star forming galaxies ( [Gallego et al. 1995 ) - 
(dotted lines). 

of measurement. Differences of at least a factor of 2 exist between different estimates of 
SFR due to these parameters. To overcome this problem, we define the luminosity density 
distribution as 

C k (P k ) = P k $k(Pk), 

where k — 1,2 correspond respectively to 1.4 GHz and Ha wavelengths. This distribution 
provides an intermediate step between the luminosity function and the star formation rate 
density and can be determined without the uncertainties due to calibration of star formation 
rate versus the luminosity. 

Luminosity density distributions for the radio and Ha samples are presented in Figures 
7(a) and 7(b) respectively. The integral of these luminosity density distributions give the 
luminosity densities, for the radio and Ha samples, as listed in Table 2. The 1.4 GHz 
luminosity density for the star-forming population in the range z < 1 (2.63 x 10 19 W Hz -1 
Mpc~ 3 ) is in excellent agreement with the value of 2 x 10 19 W Hz -1 Mpc~ 3 determined 



by Condon (|1989[ ). Also, the Ha luminosity density for the faint radio sources in this study 
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Figure 7. The luminosity density distributions, <E>(Pj.), for the 1.4 GHz radio (a) and Ha (b) surveys of faint radio sources. 
The error bars correspond to Poissonian uncertainties. 

(0.91 x 10 32 W Mpc -3 ) is consistent with the value of C Ha = 1-26 x 10 32 W Mpc~ 3 determined 
by Gallego et al. ( |1995|) , using an Ha selected low-dispersion objective-prism survey of 
nearby (z < 0.045) galaxies. Considering that these surveys are entirely independent and 
have different selection criteria, the agreement is remarkably good. 

The relative contribution from the star-forming and elliptical populations to the radio 
luminosity densities are shown as a function of their radio luminosity in Fig. 8. The star- 
forming population dominates the radio luminosity density at P1.4 ~ 10 22 — 10 23 W Hz -1 
while the ellipticals (possibly with AGN activity) contribute significantly at P L4 > 10 23 W 
Hz -1 . The implication of this result is that the main contribution to the radio luminos- 
ity density at high redshifts, using flux-limitted radio surveys, is from the elliptical/ AGN 
population (because these dominate the bright-end of the radio LF)-(Fig. 8). 

The luminosity density distribution, Ck(Pk)> can be converted to the star formation rate 
density, SFR fc , using model-dependent conversion factors corresponding to 
£1.4 



SFRx, 



8.07 x lO^WHz- 1 



Moyr- 1 
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Table 3. Integrated luminosity densities and star formation rate densities for the faint radio population. The SFR densities 
are derived from both the 1.4 GHz and the Ho measurements. The results are presented for different types of objects in the 
sample. 
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after Condon ( g99g) and 

after Gallego et al. ( |1995| ). The conversion factors and hence the resulting SFR densities 
have been adjusted to an initial mass function ip(M) oc M -2 ' 35 for 0.1 < M/M Q < 100. 
Intercomparison between SFR Q and SFR! 4 relies on the validity of the calibration of the 
1.4 GHz luminosity of galaxies in terms of their past supernova rate QCondon fc Yin 1990[ ). 



As discussed by Cram et al. (|1998|) and Cram (|1998|) , this calibration is quite uncertain. 



Using the above calibrations, the global SFR densities are estimated and listed in Table 2. 
The values of SFR1.4 are derived from the radio luminosity densities of galaxies with optical 
spectral morphologies classified as star forming, while the values of SFR Q are based on the 
Ha luminosity density for all galaxies with Ha detected. 

We find SFR1.4 = 0.032 M yr -1 Mpc" 3 for z < 1, is close to the value estimated 
from optical studies (Madau et al. 1996; Hogg et al. 1998). Also, it is similar to the value 
SFR Q = 0.031 M yr -1 Mpc" 3 we find for galaxies in the same redshift range. There is a 
rather large difference between SFR1.4 and SFR Q for the range 0.3 < z < 1. This may be 
due to a combination of several effects, including small-number statistics, greater extinction 
at higher redshift and errors arising from the use of [OII]3727 as a proxy for Ha. 
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Figure 8. The 1.4 GHz luminosity density distributions, Pi. 4 <!>(Pi.4), for the star forming population (solid line) and elliptical 
galaxies (dotted line). 

The value of SFR a derived from our estimate of the Ha luminosity density can be com- 
pared directly with that of Gallego et al. (1995) provided that we ensure that the conversion 
factors refer to the same initial mass function. Our value, SFR a = 0.032 M yr" 1 Mpc~ 3 for 
z < 1 is close to that value of 0.043 M yr _1 Mpc -3 found by Gallego et al. (1995) when 
adjusted to the same IMF. Our Ha sample is dominated by star-forming (disk) galaxies, 
implying that the star-forming radio sources are similar to those selected at optical (i.e. Ha) 
wavelengths. 

It is worth noting that the Gallego et al. (1995) sample only covers galaxies in the range 
z < 0.045. Considering this and the results in Table 2, we have not been able to establish 
evidence of cosmic evolution in the global star formation rate (Table 2), despite the fact 
that some of the radio-selected star-forming galaxies lie at redshifts z > 0.5 where the 
global SFR might exceed the local rate by a factor of the order of 10 5 ( |Madau et al. 1996| ) . 
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This conclusion reflects the statistical uncertainty of our small sample size and is consistent 
with other similar investigations (Tresse and Maddox 1998). 

4 CONCLUSION 

A sample of over 1000 faint radio sources, selected at 1.4 GHz and 'complete' (to 50%) to a 
flux density limit of ~ 0.2 mJy, is used for this study. Complete i?-band optical photometry 
is carried out to R = 22.5, with the optical spectroscopy successfully attempted for over 
225 galaxies selected at random from the radio sources with optical counterparts brighter 
then R = 21.5. The spectra have been used to derive redshifts and to make spectroscopic 
classifications of the optical galaxies into elliptical, star-forming and Seyfert classes. The 
statistical analysis is carried out on the 'complete' radio survey and the faint radio sources 
showing Ha emissions in their spectra. Conclusions from this study are summarised as 
follows: 

(i) Taking into account the various selection effects which define the two-wavelength sam- 
ple, we have derived the radio and Ha luminosity functions for both starforming and elliptical 
populations in the interval < z < 1. The LFs are well fitted to models, assuming no evo- 
lution or luminosity evolution in the form P\a{z) = -Pi.4(0)(l + z) 3 . Therefore, the faint 
radio sources are not affected by luminosity evolution at the median redshift (z < 0.3) of 
the present sample. 

Both the radio and Ha LFs agree closely with their counterparts obtained from studies 
of nearby (z < 0.05) galaxies. 

(ii) The radio (1.4 GHz) LF of the star-forming population shows a steeper faint-end slope 
than its counterpart for the ellipticals. This implies that there is a large population of faint 
radio sources undergoing star formation at z ~ 0.3, corresponding to the median redshift of 
the present sample. 

(iii) The luminosity densities are estimated and used to derive the star formation rate 
densities for both the radio and Ha data. These give a similar SFR density of 0.032 M Q 
yr _1 Mpc~ 3 at z ~ 0.3. Moreover, The SFR density estimated here is close to those found 
from optical studies. This indicates the efficiency of radio observations in selecting the star- 
forming population, free from biases induced by dust. 

(iv) The SFR density based on the Ha data, due to the radio sources in the range < z < 
1, is 0.032 M yr" 1 Mpc" 3 which is close to the local value of 0.042 M yr" 1 Mpc~ 3 derived 
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by Gallego et al. (1995). Assuming no luminosity evolution in the faint radio population, 
this suggests that the star-forming behaviour of faint radio sources is similar to that of 
star-forming galaxies selected at optical wavelengths. 

The Australia Telescope is funded by the Commonwealth of Australia for operation as a 
National Facility managed by CSIRO. LC and AH acknowledge financial support from the 
Australian Research Council and the Science Foundation for Physics within the University 
of Sydney. We thank Drs Carole Jackson and Elaine Sadler for several helpful discussions. 
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